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Abstract 

The paper is concerned with the adaptive finite element solution of linear elliptic differential equations 
using equidistributing meshes. A strategy is developed for defining this type of mesh based on residual- 
based a posteriori error estimates and rigorously analyzing the convergence of a linear finite element 
approximation using them. The existence and computation of equidistributing meshes and the continuous 
dependence of the finite element approximation on mesh are also studied. Numerical results are given to 
verify the theoretical findings. 
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1 Introduction 

We are concerned with the convergence of the linear finite element solution of elliptic differential equations 
using equidistributing meshes. An equidistributing mesh of N elements for f2 = (0, 1) is a mesh a;o = < 
a;i<---<XAr = l satisfying the so-called equidistribution principle [T2l [T9| 

J p{x)dx=^j p{x)dx, i = l,...,N (1) 

where p — p{x) is a user-prescribed, strictly positive function. Function p{x), referred to as an adaptation 
function, can be interpreted as an "error" density function, with p{x)dx being the total "error". Equation 
([ij implies that p{x) is evenly distributed among the mesh elements. 

Equidistributing meshes are known to produce optimal error bounds and have been widely used for adap- 
tive numerical solution of differential equations. Their theoretical studies have also attracted considerable 
attention from researchers; e.g., see [T2j [19l EQ] |22l [30j |38] on best approximations with variable nodes, 
[5^ HCTl m] on regression problems in statistics, [2J [MJ |33] on adaptive numerical solution of differential 
equations, and [1 H [HI EHl US HSl 123 UHl 121 [Ml ISZ] for more recent works. A focus of these studies has 
been on error analysis, i.e., to understand how accurate an approximation or a numerical solution can be on 
an equidistributing mesh. Unfortunately, this has proven to be a difficult task due to the highly nonlinear 
coupling between the mesh and the solution. The analysis can be significantly simplified by taking a priori 
meshes defined using the exact solution or some information of the exact solution. Interestingly, almost all 
of the existing analyses have been done in this way. For example, Pereyra and Sewell 34J choose a mesh to 
equidistribute a form of the truncation error and obtain an asymptotical bound for it for the finite difference 
solution of two-point boundary value problems. Qiu et al. [36l El] and Beckett and Mackenzie [H [7l [3 l29] 
investigate the uniform convergence of finite difference and finite element approximations for singularly per- 
turbed problems for meshes determined using the equidistribution principle and the singular part of the exact 
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solution. Chen and Xu |17j show that a standard finite element method and a new streamline difi^usion finite 
element method produce stable and accurate approximations for a singularly perturbed convection-diffusion 
problem provided that the mesh properly adapts to the singularity of the solution. Huang et al. |251 UHl HT] 
and Chen et al. [IS] study multi-dimensional interpolation problems using equidistributing meshes which 
depend on the function under consideration. 

The noticeable exceptions are the work [2] and f28^ where a posteriori equidistributing meshes, or equidis- 
tributing meshes determined by the computed solution, are considered. More specifically, Babuska and 
Rheinboldt [5| consider the linear finite element solution of a one-dimensional elliptic problem and develop 
a functional from a residual-based a posteriori error estimate in lieu of asymptotic approximation and coor- 
dinate transformation. Using the optimal coordinate transformation obtained by minimizing the functional, 
Babuska and Rheinboldt show that a mesh is asymptotically optimal if the residual-based error estimate 
is evenly distributed among the mesh elements. Kopteva and Stynes |28j study an upwind finite difference 
discretization of one-dimensional quasi-linear convection-diffusion problems without turning points and de- 
velop a convergence analysis for the discretization where the mesh is determined by the computed solution 
through the equidistribution principle and the arc-length adaptation function. 

In this paper we are concerned with convergence analysis for the finite element solution using a posteriori 
equidistributing meshes. The goal is to develop a systematic approach for defining these meshes such that 
both their error analysis and computation can be done in an a posteriori manner. At the same time we 
would like the approach to be general enough so that it can apply to any standard finite element method 
and have no essential limitations for multi-dimensional generalizations. Furthermore, the approach should 
be mathematically rigorous. Particularly, it should not rely on asymptotic approximation or continuous 
coordinate transformations as in [21 [57]. Several other issues, such as the existence and computation of 
equidistributing meshes and the continuous dependence of the linear finite element solution on mesh, are 
also studied in the paper. The main results are given in ^j2] 

Since Dorfler's seminal work |21j significant progress has been made on the convergence analysis of 
adaptive finite element methods based on a posteriori error estimates; e.g. see [H [T31 [HJ [Tni|3TJ |3U . It 
should, however, be pointed out that there are essential differences between those works and the current one. 
The former ones are dealt with adaptive mesh refinement using specially designed marking strategies and 
their convergence results are typically measured in terms of refinement levels, whereas the current work is 
concerned with equidistributing meshes (including their existence, generation, optimality, and error analysis) 



and our results are measured in terms of the number of mesh elements (cf. Theorems 2.1 and 2.2). It does 
not seem that the existing convergence analysis for mesh refinement can apply directly to equidistributing 
meshes and neither can the current results be covered by the existing ones. On the other hand, adaptive mesh 
refinement and equidistribution do share some common ground. For example, an equidistributing mesh can 
be generated through mesh refinement (e.g., see [lOlISS]) (and other strategies (e.g. see [24] for a variational 
approach)), and the concept of mesh equidistribution is often used in mesh refinement algorithms and 
computer codes for maximizing the efficiency of computation (e.g., see |33|V Relations between convergence 
results for adaptive mesh refinement and equidistribution may thus deserve further investigations. 

The paper is organized as follows. The description of the mathematical problem and the main results 
are given in ^ The approach for defining equidistributing meshes and analyzing the corresponding finite 
element error is developed in S|3] An iterative algorithm for computing the meshes is proposed and numerical 
results are presented in fj4] The continuous dependence of the finite element solution on mesh and the 
existence of equidistributing meshes are studied in S|5] and S|6] respectively. 

2 Main results 

We consider the boundary value problem of a linear elliptic differential equation 

-{au'Y + bu' + cu^ f, inn (2) 
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u(0) = u{l) = 0, 
where a{x), b{x), c{x), and f{x) are given functions satisfying 



and 



a{x) > ctQ > 0, c(a;) — 2^'^^^ — ^'^^ ^ 



(3) 



(4) 



(5) 



for some constant uq. Here, VF^'°°(ri) denotes the Sobolev space of functions whose derivatives are in L°°{il). 
The variational form of problem ^ and ^ is to find u £ V = Ho(ri) such that 



where 



For a given mesh 



B{u,v) = if,v), yvev 
B{u, v) = I (au'v' + bu'v + cuv)dx, (/, v) — I fvdx. 



TTh ■ a;o = 0<a:i<---<a:;jv = l 



(6) 
(7) 
(8) 



with h = max.i{xi — Xi-i), a linear finite element approximation to the solution of ^ is defined as Uh e Vh 
satisfying 

Biuh,Vh)^if.Vh), yvh^Vh (9) 



where the linear finite clement space is given by 

Vh = span{0i, (t>N-i}, 



(10) 



with (f>i^s being linear basis functions associated with mesh points x^'s. 

We are concerned with adaptive finite element solution of ^ and ([s]) using equidistributing meshes. For 
this purpose, we choose the mesh according to the equidistribution principle ([T]) or 



p{x)dx 



N' 



1 = 1, ...,N 



where 



i 



i = 1 N 



ah 



-\ 3 



rh = f + a'u'h - bu'f^ - cuh, 



{rh 



\rh\ dx 



(11) 

(12) 
(13) 

(14) 

(15) 
(16) 

(17) 



Here, Uh is the solution of ([9]), Vh is the residual function, p{x) is the piecewise constant adaptation function, 
and {rh)i is the average of rh over (x^-i, Xi). 

The choice ( 14 1 for the adaptation function is based on an a posteriori error estimate for the linear finite 
element solution; see |3] Clearly, the choice depends on the computed solution and thus the mesh is a 
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Figure 1 : Illustration of the iterative solution procedure for the finite element solution using equidistributing 
meshes. 



posteriori. Unfortunately, this also means that the mesh and the computed solution are coupled with each 
other. The system for uy^ and tt^ consists of algebraic equations (|9| and (11) and the boundary conditions 
xq = and xat = 1 and is typically solved iteratively; see Fig. [T| An algorithm of this type is given in fjl] 
The existence of the equidistributing mesh is stated in Theorem |2.4| below. 

Notation. The L^-norm on is denoted by || • ||n and other L^-norm by || • ||Lp(n), with the latter being 
extended to the situation < p < 1. Let 



(18) 



We use C as a generic constant which may have different values at difi'erent appearances. In most part of 
this paper, constants are considered as numbers that may depend on the domain and coefficients a{x), b{x), 
and c{x) of differential equation ^ but not on the solution u, the right-hand sider /, and the mesh employed 
in the finite element solution. The exceptions are Theorems 2.3 and 2.4 and f|5]and Sj6] where constants may 
further depend on u and /. 
Define 



Po = 



1 + 7 



L3 (O) , 



(19) 



where 7 is a positive constant dependent on the domain and coefficients of equation ([2|. The definition 
of 7 is given in the proof of Lemma 6.1 The same lemma also shows that po is the upper bound on the 



adaptation function defined in ( 14 1 



From ([8|, the mesh tt/j corresponds one-to-one to the {N — l)-component vector X = (xi, a;jv_i)^. 
For this reason, occasionally X is directly referred to as a mesh. Let 



'AT 



{ Xe^ 
1 



Xo 

2 

< hi< 

poN - iV' 



< xn~i < 



1, 



It is easy to verify that Sn is a closed, convex subset of 3? 
viz., 

Il^lloo = maxixil, WXeSN 



< xi < 

i-l,...,iv|. (20) 
, The set is equipped with the maximum norm, 



(21) 



This set plays an important role in the study of the existence of equidistributing meshes and the continuous 
dependence of the finite element solution on mesh. Any A^-cell equidistributing mesh is a member of this 
set. 



The main results of this paper are summarized in the following four theorems. 



Theorem 2.1 (Convergence for equidistributing meshes) Define p and as in (14-) and (15), respec 



tively. For any equidistributing mesh satisfying {11), the error for the linear finite element solution |£p is 
bounded by 



Uu^UhYWn < 



N 



(22) 
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where has the property 



lim ^foh = 1 1 r 1 1 2 



and r is the continuous "residual" function r defined as 

r ^ f + a'u' — bu' — cu. 
If further r satisfies r' E L^{il), then there exists a positive constant c such that for N > c, 



C \ 3 

N 



N 



L3 (n) 



c|k'IUi(n) ^ ' 
N 



(23) 
(24) 

(25) 



as 



The proof of this theorem is given in §3.4[ The theorem shows that the error has the asymptotic bound 



(26) 



lnn^N\\{u-u,y\\n<CM^.^^^. 



This is compared with the error bound for a uniform mesh (cf. Lemma 3.3 I 

Ci\\r\\n 



Wiu-UhYWn < 



N 



(27) 



Since 



and particularly, the left-hand side is much smaller than the right-hand side when r (~ aw") is non-smooth, 
the theorem implies that the error bound for an equidistributing mesh can be much smaller than that for a 
uniform mesh. This explains why an adaptive mesh often produces a more accurate solution than a uniform 
one when the solution is non-smooth. 

In practice, it is more realistic to use a quasi-equidistributing mesh than an exact one. A quasi- 
equidistributing mesh is a mesh satisfying 



Npihj 



< K 



z = l,...,7V 



(28) 



for some positive constant k independent of i and N . The following theorem, proved in |3.4[ shows that for 
small K, a quasi-equidistributing mesh leads to a comparable error bound as an exact equidistributing mesh. 



Theorem 2.2 (Convergence for quasi-equidistributing meshes) Define p and ah as in and (15) 



respectively. Then for any quasi- equidistributing mesh satisfying (28), the error for the linear finite element 
solution ^ is bounded by 



{u - UhYWn < 



N 



(29) 



where ah satisfies (23). 

If further r satisfies r' e L^{Q), then there exists a positive constant c such that for N > kc, 



N 



\r\\ 2 , , < A/Oft < 1 — 
"L3 (O) — V 1 _ 1 



N 



L3 (O) 



N 



(30) 



where r is defined in (24) 



Theorem 2.3 (Continuous dependence of finite element solution on mesh) Assume that f E L°°{il) 
and u E H^(fl). Then for any meshes X, X E Sn satisfying 



(31) 
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the corresponding linear finite element solutions, and uj^, satisfy 

\\{uh-u~J\\mn) < CN\\X-X\\^, 
\\uh-ujn < CN-^^\\X-X\ 
\\\uh\\1-\K\\1\ < CN-i\\X-X\ 
E denotes the energy norm associated with the bilinear form B{-, viz. 

1, 



~ i 

I 2 

I OO : 

~ i 

■| 2 

I OO : 



(32) 
(33) 
(34) 



av'^ + (c- -b')v^ 



dx. 



Theorem 2.4 (Existence of equidistributing meshes) Assume that f e L°°{^) and u G H^{il). For 
sufficiently large N (i.e., N > Nq where Nq is defined in Lemma 6.1), there exists at least an equidistributing 
mesh satisfying (11). 



The above two theorems are proven in f|5] and f|6] respectively. An iterative algorithm for computing 
equidistributing meshes is proposed in f|4] The numerical results presented in S|4] demonstrate that the 
algorithm converges for sufficiently large N and faster for larger N. This is consistent with what observed 
by Pryce |35j and Xu et al. 45J on the convergence of de Boor's algorithm for generating equidistributing 
meshes for given adaptation functions. 



3 Error analysis for finite element solution using equidistributing 
meshes 



In this section we present an error analysis for equidistributing meshes satisfying (111 and quasi-equidistributing 



meshes satisfying (28). The approach we use consists of three major steps, deriving a residual-based a pos- 



teriori error estimate, defining the adaptation function (14 1 based on the estimate, and developing the error 
analysis for the corresponding equidistributing mesh. This approach shares some similarity with that used in 
[2511261 127] for analyzing interpolation error in multi-dimensions. The main difference lies in that the current 
analysis is based on an a posteriori error estimate and is mathematically rigorous, whereas the analysis in 
[251 [57] is based on interpolation error bounds (depending on the exact solution) and valid only in an 
asymptotic sense. 

3.1 Preliminary results 

For completeness and for easy reference we list here some preliminary results without giving their proofs. 
These results can be found in most finite element textbooks, e.g. [TTl [T5] . 

Lemma 3.1 The bilinear form B{-,-) defined in has the properties 

aAv'\\l<B{v,v)<C\\v'\\l, ^veV (35) 
B{u,v)<C\\u'\\n\W\\n. ^u^v gV. (36) 

Moreover, the solution of the problem ^ satisfies 

\\u'\\n<C\\f\\n. (37) 

Lemma 3.2 Given a mesh Wh, denote by Hh the operator for piecewise linear interpolation, i.e., 

N 

U,,v{x) = V« e H^n) (38) 
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where 's are the linear basis functions associated with mesh points Xi 's. Then, for any Ki ~ Xi), 

Wv-UhvWk^ < CKWiv-UhvYWK,, yveH\K,) (39) 



\\iv~U^,vY\\K, < CK\\v"\\k^, yveH\Ki 
\\{v-Tlhv)'\\K, < C\\v'\\k„ VveH^K,). 



(40) 
(41) 



The error for the finite element solution Uh, eh — u — Uh, satisfies the orthogonality property and the 
error equation, viz.. 



B{eh,v)^{f,v)-B{uh,v), ^veV. 



(42) 
(43) 



Lemma 3.3 The finite element solution defined in satisfies 

\\<\\n < CWfWn. (44) 
Moreover, if the solution of the continuous problem |^ satisfies u G H^{il) and the mesh has the property 

C 

h = max hi < (45) 



N 



for some positive constant Ci, the error is bounded by 

\\{u-Uh)'\\n < 



CM 

N 



(46) 



Obviously, a uniform mesh satisfies the condition ( 45 1 . The finite element error for a uniform mesh can 



thus be bounded as in (46 1. 



3.2 An a posteriori error estimate 

We now derive a residual-based a posteriori error estimate for the finite element solution. The general 
procedure for this type of error estimation can be seen, e.g., in [Tl 151 HI [51 



Lemma 3.4 The error Ch = u ~ uu is bounded by 

\\{u-un)'\\l<Cr,l^CY,hnrhrK^, 



(47) 



where r^ is defined in (16), i.e., r/,, = / + a!u[^ — hu'f^ — cuh. 



Proof. Using orthogonality property ([9]), error equation (42 1, integration by parts. Lemma 3.2 and 
Schwarz' inequality, we have, for any v d V, 

B{eh,v) = B{eh,v - Uhv) 

= {f,v -U^v) ~ B{uh,v -Ufiv) 

= / rk.{v - Ilhv)dx 

< ^\\rh\\Ki\\v -Iihv\\K, 

i 

< Cj2\\rh\\Kr h^\\v\\K, 

i 



< 



Then (47) follows by taking v — in the above inequality and using Lemma 3.1 



7 



3.3 Determination of optimal adaptation function 

Up to this point the mesh has been assumed to be arbitrary and Lemma [3. 4| has been obtained for this general 
mesh. From now on we shall focus on equidistributing meshes determined according to the a posteriori error 
estimate (47 1. 

As we can see from (jlj , the key for the determination of equidistributing meshes is to define an appropriate 
adaptation function p ~ p{x) > 0. To this end, we regularize 77^ in (47 1 with a positive constant > Q (to 
be determined), i.e.. 



< 



i 

hi (an + {rh) 

i 

o^hY^hl (1 + - 



Notice that fjh r]h as a/^ ^ 0. Moreover, 



ah 



E 



h,ll+ — {rhf, 
ah 



-\ 3 



(48) 



> 



ah 

7V2 



^ V ah 



{rh 



(49) 



adaptation function p(x) is chosen as in (14 1, i.e. 



It is easy to see that the equality in the above inequality holds for any equidistributing mesh when the 

(50) 



ah 



l + — {rh)t] , i = h...,N. 



The bound in (49 1 is not the lowest due to its mesh dependence. Nevertheless, we may expect 

n 3 



ah 



a 



, 1 . 
1 + -r' 
n \ a 



dx 



as N —!■ oo. 



where a = limN^oc ah- When this is the case, (49 1 is an asymptotically lowest bound for fjf^, and in this 
sense the choice ( [SO] ) is asymptotically optimal. 

From ( 11 ), ( |49[ ), and Lemma 3.4 it is easy to see that the error on a mesh equidistributing the so-defined 
p is bounded by 



\\{u-uhy\\i<cfh = 



Cohol 



(51) 



To complete the definition, we need to determine the parameter ah- We follow [53] to choose it such that 

(Jh = Y^ < 2. (52) 

i 

In this way, roughly fifty percents of the mesh points are placed in the region where p ^ 1 |23j . From 
Jensen's inequality and (50 1, 

.h = T.h.{l^^^irh)^' 
< Y.'^^i^ + '^h' (rh)! 

i 

= l + a-^^/i,(r,)J. 
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Then (52) holds when is chosen as m (151, i.e., 



{rh)l 



(53) 



Thus, when the adaptation function p(x) and the intensity parameter ah are chosen as in (50) and (53), 



respectively, from (51 1 and (52 1 we see that the finite element error for a mesh equidistributing p is bounded 
by 



\\{u-uuy\\i<cf,i< 



Cat 



(54) 



The boundedness of as — > oo is investigated in the next subsection. 

3.4 Convergence for equidistributing and quasi-equidistributing meshes 



We notice that the adaptation function defined in (50 1 satisfies pi > 1, « = 1, A^. As a consequence, (52) 



imphes that the equidistributing mesh has the property (45) with Ci = 2. Combining this with Lemma 
we have the following theorem. 



3.3 



Theorem 3.1 Assume that p and are defined as in (50) and (53), respectively. If u ^ H (CI), then 



for any mesh equidistributing p the error in the finite element solution Uh to problem ^ is bounded by (46) 

C\\r\\n 



i.e 



{u - Uh)'\\n < 



N 



(55) 



As mentioned before, the error bound for a uniform mesh also has the same form given by (55). Although 



a bound like ( 55 1 for an equidistributing mesh is useful in some situations such as in proving the existence 



of equidistributing meshes in Lemma |3.8[ it does not show any advantage of using an adaptive mesh over a 



uniform one. In the following we shall derive a sharper bound based on the a posteriori error bound (54) 
The key is to estimate ah , and that is done in a series of lemmas. 

Lemma 3.5 (Power Inequalities) 
(i) Given a real number < q < I, for any x, y £ 3?, 



\y\' 



||a;|9_|y|9| < \x~y\'i. 



(56) 
(57) 



(ii) Given a real number < q < 1, for any two functions v and w in a function space equipped with a norm 



(58) 
(59) 



Proof. From the triangle and Jensen's inequalities, (56) follows from 



Inequality ( 57 1 is obtained by combining the inequalities 



1:^1" < |(a:-2/)+yr < 
l2/r<|(y-2:)+xr<|x-yr+|xr. 



Inequalities ( 58 1 and ( 59 1 can be proved similarly. 



9 



Lemma 3.6 For any real number < q < I and any mesh Hh for U,, 



Proof. The estimates follow from 



Lemma 3.7 For any real number < g < i and any mesh tt^ /or fi, 



|2g 



-2h'^i\n\ 



l-2g||,/||2g 



Li(0)' 



where h = max,- hi 



Proof. The left inequality is a consequence of Lemma 3.6 



To prove the right inequality, define the element-wise average of v as 

1 



vdx. 



H J Ki 



Then, from Lemma 13.51 



< 



< 



< 



< 



i 
i 
i 
i 
i 



i i 
^-VK^K+Y] I \VK^-V\^^dx 



2j2hl-''\\v-VKAZ 



From the assumption v' G L^{fl), we have 



Ki 



v{x) 



v{t)dt 



H J K. 



dx 



1 



Ki 



{v{x) - v{t))dt 
v'{s)dsdt 



dx 



Ki Jt 



dx 
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Combining ( 63 1 with ((62|) and using Holder's inequality we get 



T.^-'\\v\?l-M"L^Hn) 



< 



-2911 /||2g 



< 2h^''\n\ 



l-2g|i,/||2g 



Li(n)' 



which gives the right inequality of (61 1 



Proof of Theorem 2.1, The bound (22 1 is given by (54 1 



For (251, from (22 1 and Lemmas 3.2 and 3.6 and we have 



(64) 
(65) 



Then from (53), (651, and Lemmas 3.5 and |3.6| it follows that 

1 2 2 



> Y.hh\r\\h-Y.^h\r-rH\\l^ 



> r 



L3(n) 



which leads to the left inequality of ( |25[ ) (with c — C 3 ) . From ( 45 1 , ( 65 1 , Lemmas 

i 

i i 



3.5 



and 



3.7 



we have 



5^,3 

h ' 



which yields the right inequality of ( 25 1 



We now prove (23). In this situation, r e L (51). From the above derivation we can see that, for > , 



l + CiV-3 ||r||% <al<{l-CN- 



(66) 



Since functions having derivatives in L^{fl) are dense in L^{fl), given any e > there exists a function f 
such that 

r'eL^n) and \\r - r\\n < e. (67) 
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3dx 



Then, from Lemmas |3.5[ |3.6 [ and 3.7|) we have 

i 

i 
i 



\^dx- 



1 3 - 



Irl 3 dx 



< 



in 

If ~ r\T>dx 



\f\T^dx ]+^h. 



r\\l. 



r\ 3 dx 



< CN-^\\f'\\l^,,^+2\\r-f\ 

< CN-^\r'\\l^,,^+2ei. 



Inserting this into ( 66 ) gives 



L3(n) 



< a| < ( 1 - CN-t 



2e9 



Taking hmit as iV ^ oo in the above inequahty yields 



rll ^ 2 < hm a? < \\r\\ ^ 2 

L^{n) N^oo " L3(o) 



2e3. 



Finally, taking hmit as e ^ in the above inequahty gives ( 23 1 . 



(68) 



□ 



It is remarked that we can also use the a priori error bound (46l to estimate \\{u — u/i)'||si in (64l. For 
convenience, we list the result in the following lemma without giving the detail of the proof. 

Lemma 3.8 Assume that the solution to problem satisfies u e H^{il) and the mesh has the property 
(4-5). If r is only integrable, ah has the property If further r' G L^(ri), then ah defined in {53) is 

bounded by 



L3{n) 



CM 

N 



< o^l < \\r\ 



L3(n) 



C 



+ [j^\\r'\\LHn)] + 



CM 

N 



(69) 



Proof of Theorem 2.2 When the adaptation function p and the intensity parameter ah are chosen as 

(70) 



in ( 50 ) and ( 53 1 , a quasi-equidistributing mesh satisfying ( 28 ) has the property 

2k 



h, < 



N' 



i=l,...,N. 



Moreover, from (28 1, (48 1, and (52 1 we have 



ah 



^ / KUh ^ 3 

i 

ah{nahY 



< 



8ahK^ 



The remaining of the theorem can be proven similarly as for Theorem |2.1 
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4 An iterative algorithm for computing equidistributing meshes 
and numerical examples 

We start with describing an iterative algorithm for computing equidistributing meshes. Recall that the finite 
element equation ^ , the equidistribution relation (11), and the boundary conditions Xq = and = 1 
form a nonlinear algebraic system for the physical solution and the mesh tt/j. This system is typically 
solved iteratively; see Fig. [T] An algorithm of this type is given in the following. Starting from an initial 
mesh 7rjj'\ it produces a sequence of meshes and solutions, {ttII'\ u\^'' } . 

Algorithm for computing equidistributing meshes. Given an integer > and an initial mesh 7rj^°\ 
for fc = 0, 1, ... do 

(i) Solution of the boundary value problem using mesh Trj^'^^ This step is to find u^f^^ G vl^'^^ such that 

B{ui'\vh)^{f,vh), Vi;, (71) 

(ii) Mesh generation. This step is to compute the new equidistributing mesh using the equidistribution 



relation (111, i.e., 



where 



pW{x)dx=^ai'\ z = l,...,7V-l 



(72) 



N 



Wx e K, 



(fc) 



z = l,...,Ar 



(73) 
(74) 
(75) 



Note that the left-hand side of (72 1 is a monotone, piecewise linear function of xl and an explicit 
formula can be found as 



(k) 



\N"h 1^1=1 'h Pi 



where j is the index satisfying 



< 



1=1 



N 



Moreover, Steps (i) and (ii) define a map 



(76) 



where X^^'' and X^'^'^^'' are the (A^ — l)-component vectors corresponding to the meshes tt^'^'' and T^'h^^K 
respectively; see (20). It is not difficult to see that a fixed point of this map satisfies ([ij and thus is an 
equidistributing mesh. Furthermore, the computation can be stopped when 



\K 



(fc+i) 



= max X 



(fc+i) 



x) < e 



(77) 
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or 

AT 

where e > is a prescribed tolerance, k is a number chosen to be close to and greater than one, and is 



(k) 

It is not difficult to see that Q 'i has the properties 



the so-called quality measure of equidistribution [2S] . The second stopping criterion needs some explanation. 

s the properties 

hT.Qe% = ^^ maxg(5,>l. (79) 



N 



In addition, max^ = 1 if and only if the mesh is an equidistributing mesh satisfying Thus, if the 

(k) .... . (k) 

mesh sequence ttJ^ converges to an equidistributing mesh we will have max^ Q],y\ — s- 1 as fc — > cx); and vice 



versa. This implies that (78 1 is an effective stopping criterion. 



It is interesting to point out that max^ Qg^^ actually measures how closely the equidistribution relation 



(111 is satisfied by the mesh; see [25] for detailed discussion. Moreover, by the definition (28) one can see 



that any mesh satisfying (78 1 is a quasi-equidistributing mesh. Finally, from (78 1 and (79 1 we have 



-N{K-l) + K<Qil\<K, t^l,...,N 
where —N{k— 1) + k > when k is sufficiently close to one. 

We now present some numerical results to demonstrate the convergence of the algorithm. 
Example 3.1. This example is a reaction-diffusion equation 

-eu" + u= -2e-x{l-x)-l (80) 
subject to the boundary condition (|3|. The exact solution is given by 

u = (e^^ -e^^ + e^^ - e~^) - x{l - x) -1. (81) 

1 -e ^ ^ / 

It exhibits boundary layers at both ends of interval [0, 1] when e is small. The parameter is taken as e = 10~^. 

A typical adaptive mesh and the corresponding computed solution are shown in Fig. ^a). In Fig. ^b), 
||7rjj — TTf^ llooi niaxi((5e^ j — 1), and \\{ufi — M)'||n £^re plotted as functions of the number of iteration. It 



can be seen that both (77) and (78) are effective stopping criteria and IIttJ^"^^^ — vrl^'^'Hoo and ma,Xi{Qeq^i — 1) 
converge in a similar manner. Moreover, the solution error quickly reaches its lowest level (in one or two 
iterations for the current case). 

The number of iterations required to reach the stopping criterion max^ Qe^.j < 1.01 and the solution 
error and the modified a posteriori estimator on the final mesh of each run are hsted in Table [l] The 
results show that the underlying iterative algorithm may fail for small N but is convergent for sufficiently 
large N. Moreover, the algorithm converges faster for larger TV. These results are consistent with the 
observations made in Pryce [35^ and Xu et al. |45J for the convergence of de Boor's algorithm for generating 
equidistributing meshes for a given analytical function. It can also be seen that || (u — Uft)'||si is smaller than 
the error estimator ij and both ||(u — u/i)'||f2 and fj converge in the same order 0{l/N) as iV — > oo. These 
results conform the theoretical predictions in Theorem |2 . 1 1 and |2. 2[ 



2 ' 4 V 4 
where e = 2 x lO^'^. The exact solution is given by 



Example 3.2. Our second example is a convection-dominated differential equation 

_ ,u" + (1 - Uu' +Ul-\e]u^ e-t , (82) 



(83) 
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Table 1: Example 3.1. Iter is the number of iterations required to reach the stopping criterion max^ Qeq.i < 
1.01 or the maximum allowed number (1000 is used in the computation). \\{ufi — w)'||f2 is the error obtained 
for the final mesh of each computation. 



N 


21 


41 


81 


161 


321 


641 


Iter 


1000 


1000 


39 


4 


3 


2 


\\{uh - uYWn 


3.07 


1.41 


6.86e-l 


3.39e-l 


1.69e-l 


8.46e-2 


V 


8.10 


3.65 


1.72 


8.35e-l 


4.15e-l 


2.07e-l 




0.2 0.4 0.6 0.8 1 



100 
10 
1 

0.1 
0.01 
0.001 





llpi(k+1)-pi(k)ll 




max(Qeq(i)-1) — 




llelll • 



















2 4 6 

number of iteration (k) 
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Figure 2: Example 3.1. (a) An adaptive mesh of N — 161 points is plotted on the curve of the computed 
solution, (b) The difference between consecutive meshes (Utt,^/^^^'' — Trj^'"'' ||oo), the equidistribution quality 
measure (maxi((5eg,i ~ 1)), and the solution error — are plotted against the number of iteration, 

k. 



which has the boundary layer at x = 1 when e is small. The numerical results are showed in Fig. [3] and 
Table [2] These results confirm the observations made from the previous example. Particularly, the algorithm 
converges for sufficiently large N and faster for larger N. Moreover, the semi-norm of the error converges 
in the first order 0{1/N) as N oo. 



Table 2: Example 3.2. Iter is the number of iterations required to reach the stopping criterion max^ Qeq.i < 
1.01 or the maximum allowed number (1000 is used in the computation). ||(u;,. — 'i*)'||o is the error obtained 
for the final mesh of each computation. 



N 


21 


41 


81 


161 


321 


641 


Iter 


1000 


327 


83 


9 


5 


3 


\\{uh - uYWn 


1.20 


5.07e-l 


2.54e-l 


1.20e-l 


5.95e-2 


2.96e-2 


V 


8.08 


1.30 


6.44e-l 


2.97e-l 


1.46e-l 


7.26e-2 



Example 3.3. This example has been used by Babuska and Rheinboldt [2 . It takes the form 

- {{x + a)Pu'Y + {x + ayu= f, (84) 

where / is chosen such that the exact solution of the boundary value problem (with boundary condition (|3])) 
is 

u^{x + ay - (^'■(1 - x) + (1 + aYx) . (85) 
15 




0.2 0.4 0.6 0.8 1 50 100 150 200 

X number of iteration (k) 



Figure 3: Example 3.2. (a) An adaptive mesh of = 161 points is plotted on the curve of the computed 
solution, (b) The difference between consecutive meshes dlTrj^'^^"^' — t^I^^Woo), the equidistribution quahty 
measure (maxi((5eg,i — 1)), and the solution error — are plotted against the number of iteration, 

k. 



In our computation, the parameters are taken asp = 2, q=l, r = —1, and a = 1/100. The numerical 
results are shown in Fig. |4] and Table [3j Once again, these results confirm the observations made from the 
previous examples. 

Table 3: Example 3.3. Iter is the number of iterations required to reach the stopping criterion max^ Qeq.i < 
1.01 or the maximum allowed number (1000 is used in the computation). \\{uh — is the error obtained 

for the final mesh of each computation. 



N 


21 


41 


81 


161 


321 


641 


Iter 


4 


3 


3 


2 


2 


2 


\\{uh ~ u)'||o 


2.15e2 


1.13e2 


5.73el 


2.88el 


1.44el 


7.20 


V 


3.22e3 


1.51e3 


7.41e2 


3.69e2 


1.84e2 


9.21el 



5 Continuous dependence of finite element solution on mesh 



This section is devoted to the proof of Theorem |2.3| for the continuous dependence of the linear finite element 
solution on mesh. To this end, we need to establish some error bounds in the L°° norm, which are also needed 
in the next section in obtaining the upper and lower bounds for the adaptation function. It is worth pointing 
out that in this section we do not require that the mesh be necessarily an equidistributing mesh. Instead, 



all results hold for any mesh in Sjy- Especially, Lemmas 3.3 and 3.8 are true for any mesh having property 



(451 



We shall use two different meshes, and tt^ or X and X, in this and next sections. To distinguish the 
dependence we shall denote any quantity or function (say v) associated with mesh X by v. Moreover, in 
these two sections constants are considered as numbers that may further depend on the solution u and the 
right-hand side function / (but not on the mesh). 

We start with establishing two inequalities in Lemma |5.1| and error bounds in the L°° norm in Lemma 
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0.2 0.4 0.6 0.8 



Figure 4: Example 3.3. An adaptive mesh of = 41 points is plotted on the curve of the computed solution. 



Lemma 5.1 



\v{x)\^ < / \v\\v'\dt, Vxe [ai,5i], yveH^iai,bi) 



bi 



\v{x)\^ < 2 \v\ \v'\dt + 



1 



bi - ai 



\v\^dt, yxe[ai,bi], yveH\ai,bi). 



(86) 
(87) 



Proof. These results can readily be proven using integration by parts. 



□ 



Lemma 5.2 Assume that X G Sn and u G H'^{fl). Then the finite element error can be bounded in L 
norm as 

C\\r\\n 



|e/i||L°=(n) < 
Wh\\L-^{n) < 



N ' 
VAr+TC||r||n 



(89) 



Remark. The dependence on constant po is spelled out explicitly in ( 89 1 . This is needed for the definition 
of po ; see the proof of Lemma |6.1| 

Proof. From Poincare's inequality and Lemma 3.3 we know that the error can be bounded in norm 

as 



N 



Then from Lemma 5.1 and Schwarz' inequality we get 



TV 



C\\r\\\ 
A2 



(90) 



which leads to ( 88 ) 



To prove (89), taking v ~ e'^ and (ai,6i) — Ki (i — l,...,N) in (87 1 we have 



||e',||ioo(;,^)<^||e;,||i,+2||e;,||K.h"|k. 



(91) 



Noticing that l/ipoN) < h,, we get from (90 1 and (91 1 that 



\Wh\\l-(K,) < P,N\\e',,\\l + 2\\e'Jn\\u"h 
< PoN\\e',,\\l + C\\e'^\\n\\r\\n 



< 



iPo + l)C\\r\\l 
N 
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which gives (89). 



□ 



Note that the convergence order for ||eft ||L°°(n) is not optimal in the above lemma. The optimal can be 
obtained by making use of the Nitsche trick. 



Lemma 5.3 Assume that X G Sj^ and u G H^{fl). Then the finite element error is bounded by 



\eh\\n < 



C 



le^llo < 



C 
iV' 



c c 

||e/i||L«.(o) < —J, ||e'^||Lo=(o) < -j=, 
\\uh\\L'^(n) < C, ||w/J|l=°(o) < C*: 



(92) 

(93) 
(94) 



where the generic constant C may further depend on the solution u and the right-hand side function f . 



Proof. The second inequality of (92 1 is a consequence of Lemma 3.3 The first inequality can readily be 
proven by making use of the Nitsche trick. The proof for (93 1 is similar to that for Lemma 5.2 but makes 
use of ( 92 1 . The inequalities in ( 94 1 follow from ( 93 ) , the triangle inequality, and the boundedness of u and 
u' in L°° norm. □ 

We now consider the continuous dependence of the finite element solution Uh on mesh. 

Lemma 5.4 Assume that X, X G Sn, f G L°°(fl), and u G H^(fl). Then the finite element solutions 
Uh related to mesh X and uj^ related to X satisfy 

^ 1 

\\{uh - U/O'llii + W'^h - Uh\\1i+^ T-\Uh{x^) - U^(ij) - {uh{x^-l) - M^(i,_i))P 



i=l 



N 



J2hMx,)-uj^{S:,)\^ < CN^ \\X-X\ 



(95) 



where 



Uh 



(x) = ^Uf^{S;,)(j),{x). 



Proof. We notice that the finite element solutions can be expressed as 

N-l N-1 
Uh = ^ UhiXi)(j)i{x), = ^ Uj^{xi)(f)i{x). 



i=l 



i=l 



Moreover, ([9| can be rewritten into matrix form as 

AU = F, AU^ F, 



where 



Let 



A^ (aij), A^ (dij), ttij = B{(t).j,4>i), a^j = B{(j)j,(i)^), 

U = {uh{xi), ■ ■ • ,u/i(a;jv-i))^, U = (u^(ii), • • • ,uj^{xn-i))'^ , 

F=(Fi,...,i^jv-i)^, F=iF,,---,FN-iV. F, = if, F, = (/,0,)- 



V = (Vi, • • • , Vn^i)^, Vi = Uh{xi) - Uh{xi) = Uh{xi) - uj^{xi), 

N-l 



= Uh-Uh^ ^ Vi(l)i{x) G Vh- 



(96) 
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By subtracting the second equation from the first one in (96 1, re-grouping the terms, and taking the inner 
product of the resulting equation with V , we obtain 



V^AV + V^{A ~ Ajjj = V^{F - F). 



We now estimate the terms in (97 1 separately. First, from Lemma 3.1 we have 



It is not difficult to verify that 



V^AV = B{v,v) > aoWv'W 



^ 1 



Moreover, from Poincare's inequality we have 



N-l 



i=l 
. N-l 



ih, + hi+i) h, h,+i 
^ V. + jV,.^ + — 



i=l 



Combining (|98|-(|100| we get 

V^AV 



^ aO|| ,||2 , aO|| n|2 , «0|| n|2 , Sao,, ,,,2 
> ^ll'^ llo + y Ik Wn + ^h \\n + -^11^^ \\n 

^ 1 

II , .,,,2 , II - ||2 , «0 -1-1,^ -(^ |2 , «0 



N 



(97) 
(98) 
(99) 



(100) 



Next, we estimate the term {F — F). Noticing that 4>i + (j>i+i = 1 on {xi,Xi+i) and + (j>i+i = 1 on 
(aii, ii+i), we have 



F ~ F — 



Thus, 



f4>idx - / f4>idx 



f(t>idx - / f(t>idx + / f(t>idx - / f(j)idx 



fcjjidx- / fcjjidx - / f(j)i+idx- / f(t>i+idx 



fdx 



f<t>idx 



fdx 
f(t>idxj - 
fdx - I fdx . 



f(j)i+idx- / f4>t+idx 



N 



=1 Y^i-i 



N 



V^{F~F) = J2[ f'lJrdx- f4>^dx]{V,-V,-l)+J2iV^-l-V^) fdx. (102) 
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Denote 

x~ = min{xi, ii}, xf — max{xi,ii}. 
When (xi^i^Xi) and (xi^i,Xi) overlap, we have 



f(t>idx 



f(j)idx 



■i-l 



/(jjidx + / f4>idx + / - 4)i)dx - / f4)idx - 

Xi-i J XT' J x^-. J x~ JXi — i 



jdx 



< 



|ioo(n) {\xi-i - Xi^i\ + \xi - X 
(hi - hi) 



- Xi^i) 



hihi hi 

< ll/l|L~(n) (2|a;i_i - + \xi-Xt\ + \hi - hi\j 

< 3||/||Loo(f2) - + - Xi\) . 

On the other hand, if (xi^i,Xi) and (ii^ijii) do not overlap, we have 



dx 



f4)idx - I f4>idx 



'i-i 



\f(l)i\dx+ I \f4'i\dx 



For both cases we thus have 

f(j)idx 



f4'tdx 



< ll/llL~(f2) - Xi-i\ + \xi - Xi\) . 



< 3||/||l~(o) {\xi^i - Xi-i\ + \xi - Xt\) . 



Inserting (1041 into (1021 and using Young's inequality and hi < 2/N, we get 
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We now proceed to estimate {A — A)U. We start with computing the non-zero entries of A 
Noticing that 



tVi-l 



1 



^■iVi+l - 

by direct calculation we have 

fli.i-l — 



i>i-i = 1, on {xi-i,Xi) 



+1" (Pi+iVi+i^ (Pi+i~ T^KVi+i)' , (t>i + (t>t+i = l, on{xi,Xi+i) 



[o'4'Wi-i + Hi^^'i-i + C(/)i0i_i] dx 



Xi-l 

a, 1 



h{xi 



H,i+1 



hi 2 ^ JXi-x JXi^x 

Xi^\ 



'ViC^a; + / (c - -b')(t)i4>i-idx, 



i+i X 
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Xi + l 



dx 



a, ^a^_^ r'''(c- JbO^'rfa;, 



hi h 



i+l 



(108) 



where 



Matrix A — (a^) has similar expressions. 
Using the above expressions, we have 



a{x)dx. 



V^{A-A)U 



i=l 
N-1 



Ui Qi Oi+i fli+i - Qi 0.i \T~T /flj+l ^i+i \T~T 
(t ~- + 7 ~ )Ui - (- ~-)Ui-i - (- ~ )Ui+i 



E (H^i) - b{xi))Ui-i - {b{xi+i) - b{xt+i))Ui+i 



i=l 
N-1 



N-1 



+ [(^» - a)C/.-l + {C2,^ - C2^^)U^ + {C,+ i ~ C,+i)U,+i 



h + h + h + h. 



(109) 



where 



b'(j)idx, Ci 



{c~ -b')(f>i(l)i^idx, C2,j 



^\c-h')q^^dx. 

Xi-l ^ 



Noticing that Lemma 5.3 imphes || it'- 11/^00(57-) < C, we have 



hi iti+i hi^i 



i=l 
N 



N 



i=l 
N 



' hi h. 



N N 
i=l 



From the definitions of a; and flj, 



ai{hi - hi) + hi{ai - di) 
ai[hi - hi) H = • — 



a{x)dx - 



H \JXi_i 



a{x)dx — I a{x)dx 



From this it is not difficult to obtain 



1 Ch^ 



\athi - dihi\ — < —j-{\xi - Xi\ + \x^-l - ii_i|). 

fl 
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Thus, from the fact that l/{poN) <h<hi< 2/N it follows 



i—1 



(110) 



Similarly, we have 



I/2I = 



N-l 



i=l 
^ ,1-1 



i=l 

i=l i=l * 



(111) 



The difference Bi — Bi involved in I3 can be estimated in the same manner as for (/^.' ^ f(j)dx — J^* ^ f(l>dx) 
in {Fi - Fi). We have 

\Bi - Sil < 3 F (Ixi-i - + \xi - Xi\). (112) 

Using these estimates we obtain 



N-l 



= 2 E(^^-^')(^»^-i-^i-i^i) 

i=l 
AT-l 

< - ^ F |Si - (l^il \Ui - Ui_i\ + \Vi - Vi-il \Ui\) 



i=l 
N-l 



< E '^^m' + c\\x - x\\i + ^ E /J-i^- v,.i\ 



To estimate I4, we denote 



(c- -b')(pidx. 



Like for Bi — Bi, we have the estimates 



IQ-QI < (c-6'/2)(3|x,_i-i,_i| + 3|a;, -5,1), 



Moreover, 



|CM-Ci,i| < (c-6'/2)(3|cci_i-x,_i| + 3|x, -i,|). 



C2,i = j\c-\h')cp-ldx + jy\c-\h')4>'ldx 

= J {c- -b')<Pi{l-(j)i-i)dx + J {c- -b')(t)i{l-(l)i+i)dx 







= Cl,i — 


Ci + / (c- 










= C'l.s ~ 


Ci + / (c- 


= C'o.i+l 





(113) 



(114) 
(115) 



^i+l- 
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Then, 



h = ^ [(a-a)(c/»-i-c/.)-(a+i-c,+i)(t/,-c/,+i) 



i=l 



AT-l 



N-l 



From this we get 



JV-l N 
z— 1 i— 1 * 



(116) 



Combining (109 1, (110 1, (111 I, (113), and (1161, we get 



N -I 

i—1 



N-1 



i=l 



(117) 



Finally, (95) follows from (101), (105), and (1171. 



Proof of Theorem |2.3[ From Lemma |5.4| we can see that the key to the proof of this theorem is 



to estimate \\u'j^ — w'^IIli(o) ^^nd \\uh — For this purpose, we notice from assumption (31 1 that 

\\X - X\\oQ < miii^{hi,hi} and 

(xj_i, Xi) n (xi_i, ij) 7^ 0, i=l,...,N. 

As a consequence, we can divide [xi-i,^;;] into subintervals [xi^i, x^_i], [x^__^,xY], and [a;^,a;i]. On these 
intervals ut — Uf^ expressed as 

{uh - = (^/i (^») -^i (s^) + W/i(ii-i)0j-i (z^)) " (u/j (ij-^1 )0i- 1 (a:) + M^(£i-2)^i-2(a;) 

= (uj^ixi) ~ M^(xi_i)) 0^(2:) + (u^(ii_i) - U,;(ii_2)) (t)i-i{x) 



+ {uf^{xi-i) - u~^{xi-2)) [(l>i{x) - 0j_i(a;)^ , 



= (w,-,(£,)-M^(i,_i)) (^0,(x)-0,(x)^ 
(u/« - u;,)l[:r-,,;.] = - UJ^{5:^+l)) (t>i{x) + (m^(5,_i) - "^(^i)) 4-i(a;) 



(118) 
(119) 

(120) 



Integrating \uh — uj_\^ over the subintervals and using the above expressions and Lemma 5.3 we get 

^ ■"/illLi(2:i_i,x+_i) ^ C{hi + hi^i)\Xi^i ~ Xi-i\, 

II"'' ~ "hllLi(a;+_i,:Er) < C/ii(|xi - ii| + |a;j_i - Xi_i|), 
II"'' ^ "/illLi(xr,:E,) < C(ft,i + /ij+i)|a;j - 
Summing these estimates from i = 1 to i = yields 

||^ih-^^^||Li(f2) < C\\X-X\U. (121) 
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Moreover, differentiating (118)-(120l leads to 



Integrating — u'- 1 over the subintervals and using Lemma 



5.3 



we obtain 



Wu'h-u'jLHx-.xO < C{l+^)\X^-Xi\. 

Thus, combining these estimates gives 

\\{uH-u~J\\mn)<CN\\X -X\\^. 



Inequahty (32 1 follows from the above estimate. Lemma 5.4 and the triangle inequality. 



Next, recalling from Lemma 5.4 and (121 1 that 



\\uh - UhWn <C N\\X~ X|U, \\uh - ujLHn) < C\\X - 



by Schwarz' inequality and Lemma 5.3 we have 



\\uh-u~^\\h = I \uh - Uf^l ■ \uh - uj^\dx 



< I \uh - uf^\ ■ \uh - Uh\dx + \uh - UfJ ■ \uh ~ Uf^\dx 
n Jn 

< \\uh - Uf^\\n\\uh - Uh\\n + \\uh ~ u~h\\L^(n)\\uh - u~^\\L^n) 

< CN-^ -NWX -X\\^ + CN-'^ -WX -X\\^ 

< CN-^\\X - X\\^, 



which gives ( 33 ) 



Finally, the finite element equation ^ implies that 



From ( 32 1 we have 



IhhWl ~ \\u~Jl\ - \\f\\n\\uH-uj^\\n < CN-^\\X-X\ 



~ i 

2 

CO ; 



which gives (34) 



6 Existence of equidistributing meshes 



(122) 
(123) 
(124) 



(125) 
(126) 

(127) 



(128) 



We prove Theorem |2.4| in this section. The existence of equidistributing meshes is equivalent to the existence 
of fixed points of the map Gat defined by the iterative algorithm in ^ The key is to show that Gn maps 
Sn into Sn and is continuous. 
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Recall from (72) that the mesh Y = GnX satisfies 



Vi-l 



p{x)dx^'^, i^l,...,N 



(129) 



where p{x) and ah are defined in (14 1 and (13 1 based on the solution Uh obtained on mesh X (i.e. tt/i). 

Lemma 6.1 Assume that X e Sn and u £ H^{il). Then there exists a positive integer Nq, independent 
of the finite element approximation and the mesh, such that, for any N > Nq, 



1, 



L3{n} 



L3 (O) 



Moreover, for any N > Nq, 



l<P<Po, l<<Jh<2, 



(130) 
(131) 



where Pq is a constant defined in (19). 



Proof. The existence of A^o is guaranteed by Lemma 3.8 Its independence of the finite element approx- 
imation and the mesh is clear from (69) for the situation r' E L^{il). For the situation r e L^(ri), we can 
choose e = ^Ikll^f^^-j fo'^ some value of 9 (cf. the proof of Theorem 12. 11). Then, a smoother function f which 
is independent of the approximation and the mesh can be chosen ana an inequality similar to ( 68 ) can be 
obtained. Thus, an Nq independent of the finite element approximation and the mesh also exists for the 
situation r E L^{il). 

The inequality pi > I follows immediately from the definition of p. For the upper bound of p, from ( 50 1 , 



(130), Lemma 5.2 Young's inequality, and the inequalities ||?'||l2(si) < lkllL°°(n) and TV > 1 we have 
p, = (l + a-i(r,)?)^ 



< 1 + a 



< 1 + 2 

< 1 + 2 

< 1 + 2 

< 1 + 2 



L3 (n) 



L3 (n) 



\\r 



L°°(n) 



1 



< 1 + 2|| 

|r|| 
|r|| 



L3 (n) 



L3(f2) 



L3 (n) 



\\r 



r\\L-'(n) + {b- a')lle^llioo(n) + c\\eh\\L^^n) 



N 



N 



(i + Ci + C2 + CiV^)^ 



(l + Ci + C2)^+2||r|| 



L3(0) 



L°°(n) 



1=0(0) (1 + + C2)^ + 2||r|l^|^^^||r|||^(^) (1 + Ci + ^2)^ p§ , 



where Ci and C2 denote the constants in ( 88 1 and ( 89 1 , respectively. Notice that pq >1 and thus p§ < Pq ■ 

2 ' — ' ' — ' I 
Letting 7 = 2 (1 + Ci + C2) ^ , from the definition of pg, ( 19 1, we thus have 



i + 7lkll 



L3{n) 



L°°(0) 



= Po- 

The bounds for ah follow from the bounds for p and the definitions of ah and ah- 
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Lemma 6.2 Assume that u G H^{fl). For any N > Nq where Nq is defined in Lemma 6.1 then 
GmISn) C iS/vf or G]\i '. — > Sj\j. 



Proof. Lemma 6.1 implies that for any given X £ Sn, 1 < P < Po and 1 < cr;,, < 2. From (1291 we then 
have ^ 

iVz - Vi-i) < / p{x)dx = ^ < — 



and 



iVi - Vi-i)po > I p{x)dx ='^> ^ 

J Vi-l 



N - N' 



Thus, Y = GnX e Sn for any X e Sn- 

Lemma 6.3 Assume that f e L°°{fl), u G H^{^), and X, X e Sn- Then, 



{rut - {rht\ <C(\\X-X\\^ + |1< - ^i^J|Li(K.) + W - ^^fJliHA-.)) 



(132) 



Proof. Using the notation ( 103 1 we have 



hi - hi 2 , 1 

{Th) 



hi Jxt 



(Kp-|r,~j2)dx 



1 



hi Jx- 



Xi fX 



,+ 

i-1 



1 



+ H + )\rH\'dx-^{ + 



i JX ^Xi-i 



i-1 



KVdx. 



From the assmiiption / e L°°(ri) and Lemma 5.3 we have Hr/ilj^oo^fj) < C and ll?-,- < C. It follows 
that 



< — ~hi\- ||rh.||^oo(o) + — I \rh ~ rj^l ■ \rh + rjjdx 



hi 



h, Jx+ 



+ j-{\xi - x,\ + \xi-i - x,^i\){\\rh\\l:^(n) + lk/illi-(n)) 

C ~ C 

< ^\\X - X\\oo + ~-\\rh - rjLHK,)\\rh + 
hi K 

C ~ C 

< - X\\oo + ^(IK - u~Jl^k,) + \\uh - uJlhk,)), 
iii tii 



which gives (1321. 



□ 



Lemma 6.4 Assume that f e L°°{n), u e H^i^l), and X, X e Sn satisfying \\X - X\\oo < l/{poN). 
We also assume that N > Nq where Nq is defined in Lemma\6.1\ Then, 



(133) 



/ \p-p\dx < CN\\X -X\\oo + C{N\\X - l||oo)5. 
Jo 



Proof. Using Lemmas 6.1 2.3 and 6.3 and the inequality \\uh — uj^\\L^{n) < ||(u/t — U/i)'I1li(o)i we have 



\4-al\ < J2\\h,-hi\^,^} +k\{n^)! -(rj:)f\ 

i=l 

N 

< C N\\X-X\\^+J2k\{rk)--{rfy^\ 



i=l 
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N 



\Pi - ft I 



< C N\\X-X\\^ + {J2hr\{rH)--{rfX\)i 

i=l 

< C N\\X-X\U + - ^hYh^in) + \\uh - t^fJUi(o)))^ 

< C N\\X-X\\^ + CiN\\X-X'\\^)K 

p1 + PiPt + pi 



(134) 



< 



< 



3 ah aj^ 
\ ah — a 



< 



< 



< 



3 ahaj^ 
1,1"^ -4 



"linf. + ^lir.f^-irh)'] 



3 a/jO;^ 



3ar 



~ 9 1 9 , X 2 



flkl 



L3 (O) 



C|a| - «| I + f (ll^ - ^lloo + IK - ^'hhHK.) + hh- u-Jl^k.)) ■ (135) 



It follows from (1351, Lemma 6.3 and Theorem 2.3 that 



/ \P-P\dx = ^{ + )\P- P\dx + ^{xY ~ x+_^)\pi- pi\ 

N N 
< ^li^i - ^i) + - Xt-l)]po + ^ hi\Pi - Pi\ 



i=l 



< C\al - a|| + CN\\X - X|U + C{\\{uh - u-JWmn) + \\uh - Uf^h^n)) 

< C\al~al\ + CN\\X-X\\^, (136) 

□ 



which, together with (134), gives (1331. 

Lemma 6.5 Under the assumptions of Lemma \6.4\ Gjv is a continuous map from Sn into Spf-' 

\\GnX - GnX\\^ < CN\\X - X||oo + C{N\\X - l||oo)^. (137) 



Proof. Let Y — GnX and Y = GnX. From the equidistribution relation (1291, we obtain 

rVi rVi 



pdx 











pdx ^ — (cr,, - cr,-J 



It follows from Lemmas 16.11 and 16.41 that 

\Vi -Vi\ < I / pdx\ 

< 2 / \p-p\dx 
Jo 



(138) 



< CN\\X-X\\oo + C{N\\X-X\\oo)-', 
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which gives ( 1371 



□ 



The term involving (A^||v\r — X||oo)3 in the above lemma can be dropped for a situation shown in the 
following lemma. 



Lemma 6.6 Assume that the assumptions of Lemma 6.4 hold. If further rg — minx£n\f{x)\ > and N 

(139) 



is sufficiently large, then G^r : Sn — s- is a continuous map satisfying 

\\GnX - GnX\\oo < C N\\X - 1|U 
Proof. Recall from Lemma [5.31 that 



C C 

||eft||L=°(o) < —J, ||e/J|L=°(o) < 

N 2 VjV 



There exists a sufficiently large N such that 

\rh{x)\ = |/(a;) - {b - a')u'^{x) - cuh{x)\ = \r{x) + (6 - a')eh{x) + ce/i(a;)| 
> To - (&- a')||e',J|i=o(o) - c||eh||L~(o) 



> 2^0, 



Vx e n. 



From this and Lemma 



6.3 



we can estimate |a/^ ~ "-"^^ I ^ 



N 

Eh 



N 



< C N\\X - XW^+f^h J^""^' 



4=1 



< CA^||X-X| 



where G if^, g Ki. Combining (141) with (136) gives 



\p-p\dx < C\c 







Finally, by ( 142 1 and ( 138 1 we obtain 



C N\\X - X||oo < C N\\X - X\ 



\p-p\dx<CN\\X-X\\^, 



which gives ( 139 1 



(140) 



(141) 



(142) 



□ 



Proof of Theorem |2.4[ From Lemmas |6.2| and |6.5| we see that Gn is a continuous map from Sn to Sm- 
Recall that Sn is a closed, convex set. By Brouwer's theorem, Gjv has at least a fixed point in S^. Since 
any fixed point of Gn is an equidistributing mesh, we have proven that an equidistributing mesh exists and 
is in Sn- □ 
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